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Abstract. Three-dimensional catalogues of objects at cosmological distances can potentially yield candidate topo- 
logically lensed pairs of sets of objects, which would be a sign of the global topology of the Universe. In the 
spherical case (i.e. if curvature is positive), a necessary condition, which does not exist for either null or negative 
curvature, can be used to falsify such hypotheses, without needing to loop through a list of individual spherical 
3-manifolds. This condition is that the isometry between the two sets of objects must be a root of the identity 
isometry in the covering space S^. This enables numerical falsification of topological lensing hypotheses without 
needing to assume any particular spherical 3-manifold. By embedding in euclidean 4-space, R*, this condition 
can be expressed as the requirement that M" — I for an integer n, where M is the matrix representation of the 
hypothesised topological lensing isometry and / is the identity. Moreover, this test becomes even simpler with 



the requirement that the two rotation angles. 



corresponding to the given isometry, satisfy G Z. The 



calculation of this test involves finding the two eigenplanes of the matrix M. A CNU General Public Licence 
numerical package, called eigenplane, is made available for finding the rotation angles and eigenplanes of an 
arbitrary isometry M of S"^. 
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1. Introduction 

Although interest in cosmic topology is just over a century 
i old (Schwarzschild 1900, 1998), much interest has recently 
' developed in trying to find clues to the topology of the 
. Universe, in particular due to WMAP cosmic microwave 
' background observations (Spergel et al. 2003) which have 
been found by many to have unusual statistical proper- 
ties on the largest angular scales (e.g. Spergel et al. 2003; 
Tegmark et al. 2003; Chiang et al. 2003; Copi et al. 2004), 
even though one study across a wide range of parameter 
space failed to find any significant signal (Cornish et al. 
2004). 

Moreover, both the WMAP observations and super- 
novae la redshift-magnitude tests (e.g. Perlmutter et al. 
1999) suggest that the Universe is spherical, but with a 
radius of curvature at least as large as the horizon ra- 
dius. Indeed, the analysis of Choudhury & Padmanabhan 
(2003) rejects the flat universe hypothesis in favour of a 
spherical model at either 97% or 90% confidence levels, 
depending on which data set of type la supernovae is anal- 
ysed. 

In addition, some analyses (Luminet et al. 2003; 
Roukema et al. 2004) suggest that the spatial comoving 
section of the Universe is a Poincare dodecahedral space. 



which is a 3-manifold with as its covering space, i.e. 
curvature is positive. 

This type of hypothesis implies tight constraints on 
the curvature parameters: the non-relativistic matter 
density parameter, flm, the cosmological constant (or 
quintessence constant) 57a , and their sum, iltot — ^m + 

While much interest is presently focussing on mi- 
crowave background analyses, only a relatively small error 
in the present estimates of f^tot is required in order for 
topological lensing of sub-microwave background objects, 
such as quasars, to be detectable. 

For example, in the right-hand plot of fig. 15 of 
Gausmann et al. (2001), it was shown that for the 
Poincare dodecahedral space and {^m = 0.35, fl\ = 
0.75), i.e. for iltot — 1-1, a spike in the pair separation 
histogram (PSH) would be expected for a catalogue con- 
taining objects in the 1 < z < 3 redshift range — the 
range of most interest for large quasar catalogues. 

Is 57tot = 1.1 reasonable? While observations do seem 
to be converging closer to unity than this, there are still 
many uncertainties present. 

For example, Myers et al. (2004) found evidence for an 
extended Sunyaev-Zel'dovich Effect in WMAP data out to 
about 1 degree from centres of z < 0.2 clusters of galax- 
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ies, and they note that the contribution from clusters at 
z > 0.2 could significantly contaminate the I — 220 first 
acoustic peak, so that the uncertainties stated in estimat- 
ing iltot ^ 1; e.g. Ar^tot — 0-02 (Spergel et al. 2003), may 
be considerably underestimated due to sources of system- 
atic error such as the Sunyaev Zcl'dovich Effect. 

Hence, three-dimensional methods of detecting cosmic 
topology still remain of interest, even though consistency 
with WMAP analyses would be required for any serious 
3-manifold candidates. 

Given the many systematic effects, both due to obser- 
vation and to physical aging of objects, in catalogues of 
extragalactic objects, it is quite possible that a real sig- 
nal could be present, but that the statistical properties 
of the catalogue would be insufficient to establish statis- 
tical significance from general statistical properties of the 
catalogue. 

Any additional tests which could either falsify or sup- 
port a given topological lensing hypothesis would, there- 
fore, be useful. 

In this paper, a simple property of topological lensing, 
which is only valid for the spherical case, is presented as 
such a test. 

This is the necessary condition that the mapping or 
isometry between the two sets of objects, of which one set 
is the topological image of the other, must be a whole num- 
ber root of the identity in the covering space, . This is 
because is finite, while R'^ and are not. 

The covering space, S^, M'^ or iJ^, (which can be phys- 
ically thought of as the apparent space) relates to the 3- 
manifold, M, (which can be physically thought of as the 
fundamental domain with glued faces), by M = S^/T, 
K^/r or /T respectively, where T is a group of isome- 
tries, which can be more generally called a holonomy 
group} 

The reason why the condition that the mapping must 
be a whole number root of the identity is necessary in 
but not in flat or negatively curved space can equivalently 
be understood as the fact that the holonomy group F is 
finite for , but infinite in the other two cases. 

Testing this property enables numerical falsification of 
a topological lensing hypothesis without needing to make 
any assumption on which spherical 3-manifold would cor- 
respond to the isometry which has been found. This is 
a practical advantage for observers interested in making 
simple analyses of observational data. 

Moreover, as mentioned by Gausmann et al. (2001) in 
eq. (15) of their paper, by embedding in K^, the isom- 
etry can be written as a combination of two simultaneous 
rotations in orthogonal 2-planes in M**. Hence, the rota- 
tions in both 2-planes also need to be whole number roots 
of the identity in , i.e. their angles need to be whole 
number fractions of 27r, since the covering space can only 
be covered once. 

^ See http:7/en.wikipedia. org/w iki/Holonomy) for a for- 
mal definition of holonomy group. 
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This implies a simple and rapid test for testing whether 
or not a numerically known isometry is a root of the iden- 
tity. 

To make this paper self-contained, the complete al- 
gebraic formulae needed for these calculations, given the 
celestial positions and redshifts of corresponding objects, 
are presented. 

In Sect. 12.11 a reminder on the terminology regarding 
different types of pairs of topologically lensed objects is 
presented. 

The condition that the isometry is a root of the iden- 
tity is presented algebraically and discussed in Sect. 12.21 
An abstract representation of the isometry is written in 
Sect. 12.2.11 while the four-dimensional matrix representa- 
tions of pairs — embedded in Euclidean 4-space, R** — and 
the resulting calculation of the isometry given the objects' 
sky positions and redshifts, is presented in Sect. 12. 2!^ The 
condition that the isometry is a root of the identity is then 
presented in a form easy to calculate in Sect. 12.2.31 

In Sect. 12.31 the condition that the two rotation angles 
of the isometry are whole number roots of 27r is presented. 
In Sect. 12.3.11 the ideal case in which the isometry hap- 
pens to be expressed in a convenient basis is presented. 
In Sect. 12.3.21 the more realistic case of an isometry ex- 
pressed in an arbitrary orthonormal basis is presented as 
an eigenplane problem (in analogy with eigenvectors): an 
Earth or Solar System or Galaxy based coordinate sys- 
tem is unlikely to be aligned with the eigenplanes of the 
isometry of comoving space — in the case that the co- 
moving space we live in is a multiply connected, spherical 
3-manifold. 

A method of finding the eigenplanes is presented in 
Sect. 12.3.31 and a numerical implementation of this is 
shown in Sect. 12.4.11 

A summary is presented in Sect. |3| 

For review papers on cosmic topology, see, e.g. 
Lachieze-Rey & Luminet (1995); Luminet (1998); 
Starkman (1998); Luminet & Roukema (1999); workshop 
proceedings are in Starkman (1998) and following arti- 
cles, and Blanloeil & Roukema (2000). For comparison 
and classification of different observational strategies, see 
e.g. Uzan et al. (1999a); Luminet & Roukema (1999); 
Roukema (2002); Rebougas & Gomero (2004). 

2. Definitions and calculation of the isometry 

2.1. Cosmic crystallography: local isometries (type I 
pairs) vs generator pairs (type II) 

The isometry between two images of a single region of 
physical space yields two types of pairs of objects which 
can reveal the isometry, and uncorrelated pairs: 

(I) local or type I pairs which should occur multiple times 
independently of curvature, 

(II) generator or type II pairs, which only occur multi- 
ple times for a 3-manifold in which there are Clifford 
translations, i.e. in the flat and spherical cases, 
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Fig. 1. Upper plot: Example of A*' = 2 images of a funda- 
mental domain, where object a; is mapped to object bi, i.e. 
g{ai) = bi. Each "local" (comoving spatial) geodesic between 
"local" objects should occur N = 2 times, i.e. twice in this 
example . This frequency of occurrence is higher than for a 
Poisson distribution, in which any geodesic joining two objects 
should occur only once, not A'' times . Lower plot: Between 
the same two images of a fundamental domain, the geodesic 
connecting corresponding objects is a generator g of the (phys- 
ical) 3-manifold from the (apparent) covering space. If jr is a 
Clifford translation, then it occurs many times with (ideally) 
exactly the same length, again more frequently than that of a 
Poisson distribution. 

(Ill) uncorrelated pairs, which can be called "type III" 
pairs. 

These are schematically illustrated in Fig. 
Correspondingly, principles of detecting these pairs 
have been developed: 

(I) The type I pairs can be collected, more generally, as 
local n-tuplets (Roukema 1996), or the two-point auto- 
correlation function of the two-point auto-correlation 
function, known as "collecting correlated pairs" (CCP) 
Uzan et al. (1999b, 1999a) can be used. Gomero et al. 
(2002) noted these pairs as small deformations of the 
pair histogram, used for detecting type II pairs. 

(II) In a pair separation histogram (PSH), type II pairs 
should, ideally, show up as sharp spikes — this is "cos- 
mic crystallography" (Lachieze-Rey & Luminet 1995; 
Lehoucq et al. 1996; Gomero et al. 2002). 

As in most of observational cosmology, observed cat- 
alogues of objects are never as simple as could naively 
be hoped for making cosmological measurements. It could 
realistically be the case that a catalogue of objects con- 
tains a real topological signal, e.g. as in the right-hand 
plot of fig. 15 of Gausmann et al. (2001), where {flm — 
0.35, ft\ — 0.75), i.e. fltot = l-l, detected either thanks to 
collecting together type I pairs or type II pairs, but that 
evolutionary and selection effects make the signal of only 
weak statistical significance. 

While it is certainly possible to simply ignore such 
isometries which cannot be shown to be significant, it 



would good to be able to have some test which relates to 
their immediate geometrical properties and not to their 
membership of a statistical class. 

This is the case presented in this paper: whether the 
suspected isometry is due to type I or type II pairs, if it 
is an isometry for the spherical case, then there does exist 
a simple test capable of falsifying the topological lensing 
hypothesis. 

Let us use the notation shown in Fig.^ where only one 
suspected realisation of the isometry g is shown, mapping 
four points in one image of the fundamental domain to 
four points in another image: 



g : {ai}i=i^4 {hi}i^ 



1.4- 



(1) 



Although only three points are necessary to uniquely 
define an isometry four points are necessary if we wish 
to use the embedding in four-dimensional euclidean space 
as discussed below in Sect. 12.2.21 Given empirical uncer- 
tainties, it is probably useful to have the extra information 
provided by the fourth pair, which in the absence of ob- 
servational uncertainties and numerical errors, would be 
redundant. The case of more than four points is a gener- 
alisation of this, including more redundant information. 

2.2. Root of the identity in 
2.2.1. Abstract representation 

We can then write the first condition in the case of 
which must be satisfied if g really is an isometry of the 
covering space which generates a 3-manifold, and not just 
a coincidence: 



5"(x)^.g(.g(...g(x)...))=/ 



(2) 



where / is the identity mapping, must be true for some 
whole number n e Z. 

Moreover, the value of n should not be so high that 
"goes past" one loop of 27r around the hyper-sphere and 
evaluates to the identity after tiling the whole covering 
space twice or more. In other words, the tiling of by 
copies of the fundamental domain should only cover S'^ 
once. 

To write this more formally, using a formalism which 
also enables the expression Eq. Q , in a way that is concep- 
tually simple and easy to calculate, it is useful to embed 

in four-dimensional euclidean space, K^. 



2.2.2. 4-D representations in and calculation of the 
isometry 

By embedding in four-dimensional euclidean space, W^, 
the reader's intuition of S*^ embedded in three-dimensional 
euclidean space, can be used. 

The distance between two points (object locations) in 
comoving space can then be thought of as an arc-length 
in ]R , along the 3-sur face (e.g., Roukema 2001). 
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Let us write the j-th points of the i = 1, 2 members of 
the j-th pair of corresponding objects, i.e. where 



as 



( \ 

V ) 



(3) 



(4) 



evaluating these from, e.g. eq. (33) of Roukema (2001), 
where the objects are located at celestial positions a^j, bij 
and redshifts s, 



R 



1-0.5 



c 

Ho 



da 



— R ain{xij/R) cos Sij cosa^ 
= R ain{xij/R) cosSij sinckij 
= R sin{xij/R) sin Sij 
= R cos{xij/R) ■ 



(5) 



Since we are dealing with spherical spaces here, fi^ is pos- 
itive. 

In Mf^, the isometry g between aj and hj is a rotation 
about the origin (0,0,0,0) — is the 3-sphere (hyper- 
sphere) of radius R centred on the origin, without loss of 
generality. This rotation can be written as a 4 x 4 matrix 
M of unity determinant. 

Since M must map a.j to hj for each j (Eq. , we have 



M B 

where A and B arc the matrices 

A EE (ai Si2 as 34) , B = (bi b2 bs b4) , 



(6) 



(7) 



If the four-vectors in R"* of all four points (equiva- 
lently, of b-,) are linearly independent, then A is invert- 
ible, and multiplication by A^^ from the right-hand side 
of Eq. El gives 



M = B A- 



(8) 



which should be a matrix of nearly unity determinant, 
apart from the positional uncertainties which are consid- 
ered "acceptable" for the calculation. Discussion of what 
precision is "acceptable" is postponed to Sect. l2.4.^ below. 

Hence, M can be calculated from the celestial posi- 
tions and redshifts of the two pairs of four corresponding 
objects, together with the values of the curvature param- 
eters for which the isometry was found. 

M is a matrix representation of the generator g. 

However, even if all four points are distinct, their four- 
vectors in are not necessarily linearly independent (e.g. 
3 vectors could be aligned in the same 2-plane). Moreover, 
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if the a.j are nearly, but not quite, aligned, then the system 
could be ill-conditioned, i.e. be highly sensitive to small 
errors. Any numerical application of this method should 
either ignore choices of 4-tuplets which are not linearly 
independent enough or at least flag them and warn the 
user. 

2.2.3. Euclidean 4-D representation of the root of the 
identity constraint 

The constraint on the generator 5 = M in Eq. (|2J) can now 
be rewritten 



(9) 



Clearly, naively testing this numerically to arbitrarily 
large n would not be a practical way to test this. This 
is because in the case of incorrect hypotheses, a computer 
making numerically exact (rather than approximate) cal- 
culations would require a (countably) infinite number of 
calculations in order to check that no n is large enough to 
yield M" = /. Real computers make approximations and 
are subject to numerical errors — as n gets bigger, these 
errors would increase without limit and make the calcu- 
lation meaningless at some large value of n, unless the 
algorithm recalculated each successive estimate of M" to 
higher and higher precision, at the expense of increasing 
the computing time per calculation as n gets higher, en- 
suring no possibility of a convergent series for the total 
computing time requiring. 

2.3. Two rotation angles as roots of the 2-D identity 

2.3.1. Constraint represented in a well-chosen basis 

However, as described using a four-dimensional matrix 
representation in eq. (15) of Gausmann et al. (2001), there 
exists an orthonormal basis in which M can be expressed 
in the form 



M = 



cost/ 


— sm 


s\n6 


cos 6 

















COS( 

sin c 






— sin( 

cos (h 



(10) 



where the rotation angles 9 and (f) clearly have to satisfy 

(11) 



27r 27r „ 



Writing the least common multiple as 



(12) 



we then have the smallest value n such that M" = /. If 
either ^ or ^ are not integers, then clearly M is not a 
root of the identity. 



2.3.2. Eigenplanes 



However, the orthonormal basis in corresponding to a 
given astronomical coordinate system is unlikely, in gen- 
eral, to be the basis in which M already has this form. 

The representation of M in Eq. (|10l) has four orthonor- 
mal basis vectors. Let us write these as 
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and s and M s can be used to construct a second vector in 
the same 2-plane, yielding a pair of orthonormal vectors 
(s,t): 



S, t, U, V 

so that 

Afs = cos 9 s + sin 6 t 

Mt = — sin 6 s + cos t 

Af u = cos (/> u + sin (f) v 

Af V = — sin (/) u + cos 6 v 



(13) 



(14) 



In analogy with eigenvectors and eigenvalues, we can 
call this an eigenplane problem, where: 

A// [st] = [st] Ag, Af [uv] = [uv] A^, (15) 

and 



A« = 



cos f — sm f 
sin 9 cos 9 



C0S( 

sin ( 



— sm ( 
cos (A 



(16) 



i.e. the 1-dimensional (scalar) eigenvalue A of the tradi- 
tional problem is replaced by a 2-dimensional eigenrota- 
tion Kg or A^. 

Just as there is freedom up to multiplication by a scalar 
for the eigenvectors of an eigenvector problem, there is 
freedom up to rotation by an arbitrary (non-zero) rotation 
to find the eigenplane in the eigenplane problem. If [ s t ] 
is an orthonormal basis representing one eigenplane of A/, 
with eigenrotation Ag, then [ s t ] A, where A is an arbi- 
trary 2-dimensional rotation matrix, is also an orthonor- 
mal basis for the same eigenplane and same eigenrotation: 



Af [ s t ] A = [ s t ] AKe 



(17) 



In the particular case of interest here, there should 
exist two orthogonal eigenplanes in order for the isometry 
to correspond to a holonomy transformation, i.e. for it to 
be of interest for generating a 3-manifold as a quotient 
space of the 3-sphere centred at the origin of K^. 

2.3.3. Calculating the eigenplanes 

How do we find the basis vectors, now grouped into two 
pairs, (s,t),(u, v) spanning these two 2-planes respec- 
tively? 

Since eigenvalue problems are normally solved using an 
iterative algorithm, it seems natural to develop a practical 
algorithm for solving the eigenplane problem. 

Suppose that we already have one vector in one of the 
two eigenplanes Pi or i-2- Then, without loss of generality 
(wlog), we can write this 



_ A-fs - (s • Ms) s 
~ v/l-(s-A/s7^' 

where • is the inner product on R*. 

It is clear that s and t are orthogonal, since 



(19) 



s • t = 



s • A/s — s • (s • Af s) s 

^1 - (s • Af s)2 
(s-Afs)(l-s-s) 
^1 - (s • Afs)2 



(20) 



and 1 — s • s = since s is a unit vector. 

Next, we need a vector which is not in Pi. Consider 
three (orthonormal) basis vectors, ei, e2, e^, in the known 
basis, i.e. the basis in which and hj are calculated. 
Since these are orthogonal to one another, at most two of 
them can be in Pi . Choose one of these outside of Pi , let 
us call it ei wlog. Then, 



ei — (ei • s) s — (ei • t) t 
Afu- (u- Afu) u 
^1 - (u • Afu)2 ' 



(21) 



are an orthonormal pair spanning the second plane, P2- 
We then have 9 and 6 from 



cos ^(s-Afs) 
cos^"'^(u • Afu) 



(22) 



Hence, finding a solution for (s,t),(u,v) reduces to 
finding at least one vector s in one of the two planes Pi , P2 • 

The intersection of Pi (or Pj) with is a great circle, 
centred on the origin. This must intersect somewhere with 
the 2-sphere defined by 



Si 



23 



span{ei,e2,e3} n S^ . 



(23) 



Hence, it is sufficient to search through S\2z looking 
for at least one vector which lies in Pi or P2. In fact, it 
is sufficient to search through half of this region, since 
polarity is unimportant. 

To write down the condition for a vector s to be in one 
of the two eigenplanes, note that the three vectors 



s, Afs, Af^s, 



(24) 



must all lie in the same eigenplane. This eigenplane in- 
tersects with 5*^ in a great circle, i.e. rotating twice cor- 
responds to a single rotation by twice the angle of the 
original rotation. We can write this condition as 



cos~^(s • Af^s) = 2cos"^(s • Afs) 
for a rotation angle 



(25) 



s e Pi 



(18) 



9 = cos"^(s • Afs) < 7r/2 



(26) 



6 

or 
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2-K - cos"^(s • M^s) = 2cos"^(s • Ms), 
for 

■K/2<e = C0S"^(S • Ms) < TT. 



(27) 



(28) 



Hence, an iterative search to sucessively preciser res- 
olution through Si2z to find a vector s satisfying either 
Eq. H25|l or Eq. (|27|l to the desired numerical precision 
yields one basis vector s of the eigenplanes, and Eqs. 1)19(1 . 
((21(1 yield the other three basis vectors. 

Tests using the eigenplane package fSect. 1^4.1(1 sug- 
gest that the number of steps n to reach a precision in 
radians of A6' is 



n = 2 - 4 logio M 



(29) 



i.e. for M - lO^'^-lO"^, about 14-18 iterations are suffi- 
cient. 

The rotation angles corresponding to the isometry are 
those in Eq. (j^ . 

If these angles do not satisfy Eq. ^2 then the isometry 
is not a root of the identity and the topological lensing 
hypothesis is false. 

2.4. Numerical aspects 

2.4.1. Numerical implementation 

The package eigenplane, which is a free software (GNU 
GPL licence) implentation of an iteration algorithm to 
generate the two rotation angles 9, (j) (and as a side prod- 
uct, a choice of four orthonormal vectors (s,t),(u, v) 
defining the two eigenplanes), given an input isometry 
M, has been prepared and is available for download at 
|http : / / cosmo .t orun . pl/GPLdownload/ eigen/ ^ 

The package is self-contained, apart from requiring 
fortran and C compilers and the public domain lin- 
ear algebra library bias. Installation (help in README 
and INSTALL files) is by the standard ./configure; 
make; make install sequence, including standard op- 
tions such as a non-root user installation directory via 
— pref ix=PREFlX. Following installation, help is avail- 
able with the commands eigenplane — help and info 
eigenplane. 

2.4.2. What are "acceptable" levels of errors in the 
positions? 

As was discussed, e.g. in Roukema (1996) and later papers, 
these consist, in principle, both of "measurement" errors 
— to what extent the observed values of angular positions 



^ At the time of proofchecking this article, the current version 
of the package is eigenplane-0 .2.3. Users are welcome to add 
features, correct bugs and distribute modified versions. 



ctij, 5ij and redshift Sij^ incorrectly represent the true po- 
sition if peculiar velocity is ignored — and of "movement" 
error — to what extent the objects move relative to the co- 
moving reference frame between the different epochs (red- 
shifts) at which they are observed. 

In practice, the precision on the angular positions is 
negligible compared to all other errors. The former are al- 
most always more accurate than an arcsecond, which at a 
redshift of s = 2 and (fim = 0.35, 17a = 0.75) corresponds 
to a tangential error of better than Q.Q2h~^ Mpc. 

The measurement error in the redshifts, e.g. of quasars, 
can usually be obtained to As < 0.001, but in big cata- 
logues is more often As < 0.01. Again at a redshift s = 2 
and (Sim — 0.35, 17a — 0.75), these two errors lead to dis- 
tance errors of about l/i~^ Mpc and 10/i~^ Mpc respec- 
tively, at least 50 times larger than that induced by any 
angular error. 

The movement error between observations at different 
redshifts depends on how different the redshifts are, and 
on hypotheses of structure formation within the general 
model of gravitational collapse from linear perturbations. 
The objects most likely to be observed are the brighter 
ones, which tend to lie in the heavier potential wells, which 
can be expected to move more slowly relative to the co- 
moving reference frame. 

If we estimate a maximum for the peculiar velocity 
as a mean of 400 km/s then between two substantially 
different redshifts, e.g. differing by 4 Gyr, the total dis- 
placement (relative to the comoving reference frame) is 
about l.&h~^ Mpc. Moreover, simulations suggest that 
some objects may "stream" from voids towards filaments 
and along filaments towards the deepest potential wells 
in the cosmic web of gravitationally bound structure — 
i.e. movement by up to a few megaparsecs over a big frac- 
tion of a Hubble time is realistic for objects forming some 
distance away from the biggest clusters of galaxies. 

So likely errors are in the range l-10ft.~^ Mpc, both 
from redshift measurements and from possible movement 
relative to the comoving frame. 

How much do these affect testing an isometry as a root 
of the identity? 

Again, for (fim = 0.35, IIa = 0.75), the curvature ra- 
dius is i? « 10/i~^ Gpc, so the errors in three-dimensional 
position, i.e. within the 3-surface , are at the level of 
about 10"'' to 10"'^ curvature radh, i.e. about 10~^ to 
10"'^ radians. If these lead to a similar error in the es- 
timates of the angles 9 and 0, and the magnitude of 
max(0, (j)) is a big fraction of the distance to the sur- 
face of last scattering, e.g. > hh^^ Gpc ~ 0.5 radians, 
then the fractional (relative) errors in max(^, ^) should 
be of the same order of magnitude, about 10^** to 10"'^. 
A reasonable lower bound could be placed on max(^, </>), 
for a given (f2ni, 17 a) pair, in order for the isometry to 
be realistic (e.g. given cosmic topology constraints from 



^ The usual variable z is used here for positions in ' 
write s (red- "shift" ) instead. 
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the cosmic microwave background), removing isometries 
where max(0, 0) is too small. So, for holonomy groups of 
order up to ^ 10^, the probability of false (chance) isome- 
tries yielding integer solutions is about 10~'^ to 10^^. For 
larger orders, the probability of false solutions necessarily 
becomes higher. 

It remains possible that the smaller angle, min(0, 6*), 
could be quite small and would have a larger relative er- 
ror (e.g. 10~^-10~^), so in many cases would constitute a 
weaker test of candidate isometries. 

If the values {ilm, ^a) are decreased to approach ^lm + 
fl\ = 1, then the curvature radius R increases rapidly but 
the distance to the surface of last scattering only changes 
slightly, so max(6', 0) decreases. This causes max(^,^) 
to increase rapidly, so even while the relative error will not 
increase, the absolute error will increase to above ±1, in 
which case integer solutions will necessarily be found for 
chance isometries. 

However, as pointed out by Mota et al. (2005), as the 
curvature radius becomes larger and larger, observation- 
ally detectable isometries in spherical (and hyperbolic) 3- 
manifolds become closer and closer to those for the flat 
case, in which case methods used in the flat case are likely 
to become preferable. 

3. Discussion and conclusion 

Of course, checking that an isometry is a root of the iden- 
tity is only a necessary condition for the topological lens- 
ing hypothesis to be correct, it is not sufficient. Moreover, 
even if it genuinely is a case of topological lensing, the 
isometry might, in principle, not be an isometry between 
adjacent copies of the fundamental domain — or in other 
words, the group it generates as the n-th root of the iden- 
tity might only be a sub-group of the full holonomy group 
(Weeks et al. 2003). 

Note that there is no point checking the Poincare 
dodecahedral hypothesis of Luminct et al. (2003) and 
Roukema et al. (2004) by this method, since the isome- 
tries there are already known to be roots of the identity; 
the assumption of a particular class of 3-manifolds was an 
input assumption. 

The situations where this test can be useful are those 
where a candidate isometry is obtained for the spherical 
case without having assumed any particular 3-manifold, 
but only having assumed that curvature is positive. 

How can these candidate isometries be obtained? In 
other words, what are the known systematic methods for 
exploring a catalogue of objects extending to high red- 
shifts in order to find a suitable pair of 4-tuplets of corre- 
sponding objects with a corresponding g which should be 
tested? 

In principle, the method of finding local n-tuplets of 
Roukema (1996) provides the answer: compare all "local" 
(less than a few hundred Mpc) "configurations" of ob- 
jects against all other "local clusters" , each time checking 
whether the mapping / between the two 4-tuplets is an 
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isometry g or not. In practice, optimising the algorithm 
for application to large, recent catalogues is a non-trivial 
task, since the number of n-tuplets climbs rapidly as the 
number of objects increases: the nai've approach with a 
large catalogue quickly becomes impractical even on the 
latest supercomputers, though some suggestions for short- 
cuts to the algorithm are made in Roukema (1996). 

A related approach might be to use the "collecting cor- 
related pairs" method (CCP) (Uzan et al. 1999b, 1999a; 
Gomero et al. 2002), to first obtain a list of 2-tuplets 
(pairs) "most likely to be matched type I pairs" , then try 
to combine these into 4-tuplets, test whether or not the 
mappings / are isometries g, and finally test whether any 
"best" isometry g is a root of the identity, using the al- 
gorithm presented here. The definition of the "most likely 
to be matched type I pairs" might first be some way of 
choosing the highest signal-to-noise ratios as a function of 
the curvature parameters (rim,^^A) (as is recommended 
in Uzan et al. (1999b, 1999a)), and then for each choice 
of (flrm^A), choose the bins in the pair separation his- 
togram with the highest numbers of pairs as the "most 
likely to be matched type I pairs" . 

Whereas Uzan et al. (1999b, 1999a) hope for a strong 
signal, use of the algorithm presented here could poten- 
tially enable detection of a candidate 3-manifold even if 
the signal is weak. 

Type II pairs — generator pairs — detected as high 
spikes in a pair separation histogram (PSH) (Lachieze- 
Rey & Luminet 1995; Lehoucq et al. 1996; Gomero et al. 
2002), could also be used to generate pairs of 4-tuplets 
and isometries g. 

Another caveat is that this test necessarily depends on 
the choice of curvature parameters (f2ni,^^A)- 

If an isometry g is known numerically based on, e.g. a 
pair of 4-tuplets of observed objects for a given choice of 
(rim, ^^a), then it is fairly likely that the same pair of 4- 
tuplets will be approximately isometric for nearby values 
of (17m, ^^a), especially if the 4-tuplets are local 4-tuplets 
(composed of type I pairs rather than type II pairs). So, 
the range of (f2m, ^^a) for which the mapping between 
a given pair of 4-tuplets remains an isometry needs to be 
tested, unless external constraints on (fJm, ^a) are consid- 
ered acceptable. This (unfortunately) increases the chance 
that a false (chance) isometry will yield a valid root of the 
identity. 
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